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It is shown that in adiabatic approximation for nuclei the many-component 
Coulomb system cannot be described on the basis of the grand canonical ensemble. 
Using the variational Bogolyubov's procedure for the free energy, the Hohenberg- 
Kohn theorem is proved in the canonical ensemble for inhomogeneous electron 
gas at finite temperature. The principal difference between consideration in the 
framework of quantum statistics in the canonical ensemble and quantum-mechanical 
consideration of a finite number of particles in infinite volume is established. The 
problem of universality of the density functional for describing the inhomogeneous 
electron density in a disordered nuclei field is considered. 
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I. INTRODUCTION 

In 1964 P.Hohenberg and W.Kohn [1] proved the remarkable theorem, which initiated 
the new method in statistical physics, i.e, the density functional theory (DFT). The DFT 
is widely used, first, in the condensed matter theory (see,e.g. [2-6]). As applied to the 
investigation of properties of the real matter the DFT is based on the consideration of a 
large number of electrons and nuclei, interacting by the Coulomb law (Coulomb system, - 
CS) [7]. 
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The large difference in the electron and ion masses in the CS allows the use of the 
adiabatic approximation (Born-Oppenheimer approximation). As is known, the adiabatic 
approximation is reduced to the following procedure. At first, the electron wave function 
has to be found in the field of nuclei which have some fixed positions. The obtained wave 
functions and the electron subsystem energy depend parametrically on the coordinates of 
nuclei. In this case, according to the Hohenberg-Kohn theorem the electron system energy 
is the functional of the electron density [1]. In the next stage, the electron energy is used as 
an effective many-particle indirect interaction of the nuclear subsystem. The electron energy 
increases with the average distance between nuclei [8]. Therefore, in the electron ground 
state, electrons stabilize the nuclear subsystem. The internuclear repulsion prevents the 
whole system from collapse and an equilibrium configuration can exist. The Fermi statistics 
of electrons plays a crucial role in collapse prevention [9] . 

Application of the DFT to description of thermodynamic properties of matter was ini- 
tiated by the work of Mermin [10]. In this paper, the Honenberg-Kohn theorem for the 
nonhomogeneous electron gas at finite temperature has been proved in the grand canonical 
ensemble. However, there is a certain problem in constructing the Grand canonical ensemble 
(in what follows, the grand ensemble) for the CS. In conventional consideration (sec, e.g. 
[11]) the grand ensemble is constructed as a set of canonical ensembles with different num- 
bers of particles of various types in each canonical term. Furthermore, the summation of 
these terms over all numbers of particles leads to the grand ensemble. For such a procedure 
for the Coulomb system containing electrons (index e) and nuclei (index i) in a volume V 
at temperature T in each canonical term, the quasineutrality condition should be satisfied. 



where nr =< Na >(^"") /V is the average particle density, Na is the operator of the total 
number of particles of type a with the mass iria-, average density n^, and charge ^r^e, angle 
brackets < ... mean averaging over the canonical ensemble. In this case, the numbers 

of particles are not already independent, and the main advantage of the grand canonical 
ensemble, which allows calculation of the Green's functions of noninteracting particles, fur- 
thermore, construction of the diagram technique perturbation theory in the quantum case, 
disappears. Therefore, this way for constructing the grand canonical ensemble is absolutely 
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unpromissing for calculation opportunities of the statistical theory. 

In the adiabatic approximation, the situation is even worse. In this case, using the 
grand ensemble, it is in principle impossible to satisfy the quasineutrality equation ([T]) for 
each nucleus configuration, which is the necessary equilibrium condition for the CS [7]. 
This means that the conventional DFT for inhomogeneous electron g clS clS db method for 
describing the properties of matter is possible only within the canonical ensemble (at least 
in the adiabatic approximation for nuclei). 

In the present paper we generalize finite temperatures in the canonical ensemble the 
Honenberg-Kohn theorem which was established in [1] for the ground state energy of the 
inhomogeneous electron gas. 



II. VARIATIONAL PRINCIPE FOR FREE ENERGY AND THE 
HAMILTONIAN OF THE CS 

In 1956 N.N.Bogolyubov formulated the variational principle which establishes the upper 
estimate for the free energy of the system with the Hamiltonian H = Hq + Hi 

F<Fo+< Hi >r, (2) 

Fo = -T In Sp jexp | ; < //^ >r= Sp !^Hi exp (^^^^) } • (3) 

In this case, the smallness of the Hamiltonian Hi is not supposed (see also [12]). The rigorous 
proof of relation ([3]) can be found in [13]. In alternative formulation, the variational principle 
([2]) for the free energy is presented in [14,15] for the Bardeen-Cooper-Schrieffer model in the 
superconductivity theory (see also [16]). In the second quantization representation, the CS 
Hamiltonian has the form 

H'^^ = Hee + Hec + Hcc, (4) 
Haa = -^ j ^+(r)V2^(r)rfr + j M,„(ri - r2)^+(ri)^+(r2)^.(r2)^.(ri)rfrirfr2, (5) 

Mec(ri - Y2)Ne{Y-i)Nc{v2)dY-i_dT2; Uab = • (6) 
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Here Uabij) is the Coulomb interaction potential for particles of types a and h, \E'"'"(r) and 
\E'(r) are the field operators of creation and annihilation, A'a(r) = \l'+(r)\E'(r) is the particle 
density operator for the type a, which in the coordinate representation, is written as 

Na 

iV,(r) = J]5(r-Rj, (7) 

i=l 

where Ra is the coordinate of the nucleus a. It should be noted that the total number Na 
of each particle type in the canonical ensemble is a C-number, 

Na = j Na{r)dv =< Na >^'^" . (8) 

The free energy of the quasineutral CS under consideration, containing the fixed number of 
electrons N^ and nuclei Nc in a volume V is given by 

F^^ = -Tln5p{exp(-^)}. (9) 
In the adiabatic approximation the free energy F'-'^ can be written as 



TlnSp^ |exp (^^^^^J j ; F,, = -TlnSp^ |exp (^-^^^^^J | . (10) 

The function Fee (E]) is the free energy of the electron subsystem in the external field of 
immobile nuclei. Moreover, the operator of nuclei density Nc{r) in the Hamiltonian Hec ( ITOll 
is a C-number. 



III. THE HOHENBERG-KOHN THEOREM FOR FREE ENERGY OF 
INHOMOGENEOUS ELECTRON GAS IN THE CANONICAL ENSEMBLE 

Now let us show that in the canonical ensemble, the inhomogeneous density of interacted 
electrons (for the definite electron number N^ in a fixed volume V) in a certain scalar external 
field with the potential (pextiT") 

ne{r) =< Neir) >(^'^") (11) 

singly- valued determines this field (with an accuracy to an insignificant constant). Let 
n^e\r) be the inhomogeneous electron density in an external field, which is defined by the 
scalar potential (p^exti'"^) ■ this case the free energy Fe^^ of the electron system equals 

.(1), 



(12) 
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In turn, let n'^\r) be the inhomogeneous electron density in an external field, which is 

('2) (2) 

defined by the scalar potential (plJ^{r), and the appropriate free energy Fe is given by 

F« = -ri„Sp,{exp(-^^i±l4^£l^j}. (13) 

Let us assume that the same electron density can correspond to different external potentials 
n«(r) = nf\v) = ^^(r); cp^^J,ir) ^ ^fi(r) + const. (14) 
According to the Bogolyubov variational principle (EJ, (|3]) 

F^'^ < + I (r) - v.S(r)}nf (r)rfr, (15) 

Fi^) < + I {^fj^(r) - ^2*(r)}n«(r)rfr. (16) 

Summing up these inequalities (fT5|) and (fT6!) we arrive at the contradiction 

+ <F«+Ff . (17) 

This means that the initial assumption is wrong. Therefore, one can assert that the inho- 
mogeneous density ne(r) of interacted electrons (for a fixed full number of electrons in a 
certain volume V) in the canonical ensemble in a certain scalar field v^ext(r) singly- valued 
defines this potential 

ipe,t{r) = (fextir, T, V, N,, {ne(r)}). (18) 

Moreover, in the canonical ensemble the free energy (see f|T2|) and (|T3|) ) and the average 
density ne(r) (see (ITT]) ) of the inhomogeneous electron system at temperature T is completely 
defined by the total number of electrons Ne, the volume V and the external potential V5ext(i') 

F, = F,{T, V, N„ {^extir)}); ^^(r) = ^^(r, T, V, {</^ext(r)}), (19) 

with the normalization condition 



J n,{r)dr =< N, . (20) 



The above proof shows that the free energy functional Fe{{ne{r)}) of the inhomogeneous 
electron system isthe functional of density ne(r) at fixed temperature T and volume V, 

F, = F,{T,V,N,,{ne{r)}). (21) 
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It should be stressed that in general the explicit expression of the free energy functional 
Fe{{ne{i^)}) depends on the form of the external potential (y9e^t(r), and on the form of the 
interparticle interaction Ua6(r) (see (El)). This means that the "universality" of the density 
functional is related to the thermodynamic parameters of the system. Therefore, for the fixed 
external field potential and fixed interparticle interaction, the density functional for the free 
energy is invariable for arbitrary thermodynamic parameters. The exception is the new 
phase formation, e.g., liquid-solid transition, etc. The determination of the explicit form of 
the density functional (fT8|) for the external field (/9ext(r) has the same level of complexity (see 
below) as the determination of the density functional (l2Ti) for the free energy Fe({ne(r)}). 



IV. ELECTRONS IN THE EXTERNAL FIELD OF NUCLEI: QUANTUM 
MECHANICS AND QUANTUM STATISTICS 

Let us now pay attention that the inhomogeneous average density profile ne(r) is caused 
not only by the existence of the external field V5ext(r) but also by the finiteness of the system 
volume V. To eliminate the boundary effects, one should pass to the limit V oo. In this 
case the definition of the limit 



1 f 
lim — / np(r)dr 

V^oo V J 



(22) 

V 



is essential (see ( 120|) ). 

As applied to the problem of electrons in the external field of nuclei the external potential 
V?ea;t(r) has the form 

Nc 

¥'ext(r) = ^Ma6(r-R^). (23) 

1=1 

We take into account that V is very large but finite, and the realpassage to the limit V ^ oo 
is possible only after averaging over nuclei coordinates. With this remark, the limit (122|) . 
can be considered in two physically different cases: 

Variant (A): for V — )■ oo, the number of electrons remains finite. This case corresponds 
to the consideration of atoms, molecules, etc. in the framework of the traditional quantum 
mechanics in the adiabatic approximation for nuclei (see, e.g., [17]). In this case, the limit 
( 122|) equals zero, the conceptions of temperature T and, hence, the free energy become 
irrelevant. The problem is reduced to the determination of the average energy EjP^ and 
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average density ne(r) of the ground state for a given number of electrons in the field of a 
finite number of the immobile nuclei. 

Taking into account the quasineutrality condition 



= (iV„ {R^, z,) = (iVe, n,(r, {R^), z,) , n,(r) = (r, iV„ {R^, z,) .(25) 



Relations ( 125|) respond to the conventional quantum-mechanical problem of the determi- 
nation of the ground state of electrons localized on one (atom), two (molecule) or more 
Coulomb centers in infinite space. In the most general formulation the problem comes to 
the determination of the absolute minimum of the ground state and the corresponding dis- 
tances between nuclei in the many-center problem. It should be stressed that there is only 
one length parameter in this problem, i.e., the Bohr radius Oq = h'^ / (m^.e'^) . Therefore, for 
the characteristic size Lg of the density ne(r) inhomogeneity with the accuracy to the z^. 
and Ai'c-dependent constant one can write Lg ~ oq. Taking into account (El]), this means 

Variant {B): for V — t- oo, the number of electrons (as well as the number of nuclei 
Nc) proportionally increases in such a way that the average electron density n^e^ (see ([T])) 
remains a finite non-zero value. This consideration responds to the "thermodynamic limit" 
(see, e.g., [18]) and corresponds to the use of statistical (particularly, canonical) ensembles. 
Namely, for this description within quantum statistics we use the concept of temperature 
and free energy. The passage to the thermodynamic limit allows us, instead the so-called 
"extensive" thermodynamic values, in particular, the free energy Fg, to use the "intensive" 
thermodynamic functions, in particular, the free energy fy of the system under consider- 
ation, which is related to the unit volume (see, e.g., [11,18]). Then, taking into account 
(II]),(II9])-(I23]), one can write 




(24) 



a=e,c 



the expressions for Ee and ne(r) can be written as 



Fe(T,y,iV„K(r)})= / /^(T,n(°),K(r,{R,^})})t^r 



(26) 



V 




(27) 
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Therefore, the functions fy and ne(r, {R^}) depend on the thermodynamic parameters, 
i.e., the temperature T and average electron density n^e^ in the volume V of the system 
under consideration. Considering the grand canonical ensemble instead the dependence on 
the average density ni^"* we arrive at dependence on the chemical potential of electrons fi^e\ 



dF \ Iff df^ \ 



(28) 



V 

If the external field is absent (in a homogeneous system), relations (!26|) .( p7|) takes the form 

Fe(T, V, AT,) ^ (T, n(°)) V, n^ir) ^ n^^\ (29) 

For the definite form of the functional Fe{T, V, Ng, {ne(r)}), the inhomogeneous density ne(r) 
is determined, taking into account f|T8l) . from the variational equation for the free energy 
(see, e.g., [11], [17]) 

5Fe = 0, j (ne(r) - n(°)) rfr = (30) 

at fixed values of the variables T and rLe\ 

It is easily seen that, taking into account ( ITOi) . ( fT2|) . ( 1T3|) . the variational equation 

^E, = 0, j (ne(r) - n(°)) dv = 0, (31) 

directly follows from Eq. ( l30l) . Here the average energy E^ of the system under consideration 
in the external field is {He)^^^^\ where He is the Hamiltonian of the system in the external 
field. In turn, the solution of the variational equation (!3T|) leads to the conclusion that 
the inhomogeneous density ne(r) is determined from the relation ffTTl) . This statement is 
equivalent to the known statement that the variational equation (!3T|) for the energy of 
the quantum-mechanical system leads to the Schrodinger equation for wave functions and 
corresponding energy levels [17]. If we now apply the variational equation (|30l) to the right- 
hand side of the Bogolyubov inequality (|2]) we immediately arrive at the formulation of 
the Hohenberg-Kohn theorem for the free energy, which reads that the external field is the 
single-valued functional of the inhomogeneous electron density. 

Moreover, according to the above consideration, the density functional for the ground 
state energy E^^^ (l25ll of the finite number of electrons in an external field of nuclei in 



the limit V oo (quantum mechanics) is principally different from the functional for 
free energy F^ of the electron system in the external field of nuclei in the thermodynamic 
limit (quantum statistics), even for the case of strong degeneracy (T — )• 0). This principal 
difference is conditioned by different values of the limit (l22l) . 
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V. LOCALIZED AND DELOCALIZED ELECTRON STATES 

From physical reasons [7] one can suppose that in a wide temperatures T range and 
average electron densities n^e^ the electron states in the external field of a large but finite 
number of nuclei can be divided into two groups [19,20,21]: 

- the "localized" states (index "loc"), in which electrons and nuclei exist in the form of 
bounded complexes, as "atoms" , "molecules" etc. 

- "delocalized" states (index "deloc"), in which electrons are spread over the entire system 
volume. 

It should be emphasized that such "atoms" and "molecules" are not identical to atoms 
and molecules conventionally considered in quantum mechanics. The essential difference 
is conditioned by the opportunity for each of the large (in thermodynamic limit - infinite) 
number of electrons (due to their identical nature) to occupy the electron state of the fixed 
"atom" or "molecule" . This fact is refiected in the statistical description of the system 
under consideration. Therefore, the electron states in "atoms" and "molecules" depend on 
the thermodynamic parameters, i.e., the temperature T and average electron density n^e^ 
(see,e.g., [22]). Then, according to the definition ( l20l) . the average inhomogeneous electron 
density can be represented in the form 

n,(r, {R^) = n('-)(r, {R^) + '-)(r, {R^), (32) 

where ne°^\r, {R!^}) and ni^^''"^\r, {TV^}) are the inhomogeneous densities, responding to 
the "localized" and "delocalized" states, respectively. Following this division one finds 

= iV^('°^) + N^^deioc)^ j^iioc) ^ j ^(M(r, {Ii1})dr, N^/^^""^ = j n('^"'°'=) (r, {R^})c/r(33) 

Here N^'"'^ and nI,"''"''^ are the total numbers of "localized" and "delocalized" electrons, 
respectively. In the thermodynamic limit, relation ( 133|) reads 

^{0) ^ ^{loc) ^ ^{ddoc) ^ (34) 

where n^e"'^^ = Ne"^"* j V and nt'^^"'^^ = '^^^^"'^^ j y are the average electron densities in the 
"localized" and "delocalized" states, respectively. Their ratio depends on the thermody- 
namic parameters of the system. 

It should be emphasized that the electron states are considered in the irregular field 
of nuclei [19]. In this case, the role of the eigen quantum numbers is played the energy 
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of the "localized" (with the discrete spectrum) and " delocalized" (with the continuous 
spectrum) states and their spin. In general, the coordinates of nuclei R^, can be added 
to the set of quantum numbers for the localized states. The localized states with fixed 
energies can be localized at different points of the system. Therefore, one can interpret these 
"localized" states as degenerate with respect to values {R^}. However, this interpretation 
has relative character, since the vectors {Rf} and their components are not eigen functions 
of any operator, affecting on the electron variables. Strictly speaking, the coordinates R? 
are the parameters appeared due to the use of the standard adiabatic approximation, when 
nuclei are considering as immobile particles [19]. Taking into account that the electron states 
with the discrete energy spectrum can be localized on one, two or more nuclei, the electron 
density ne°^\r, {R^}) can be written as 

^Hr, {R.}) = E ^(iT^(r - + E E ^(2TV, R^) + • • • , (35) 

where n^^^"^^^ and n^l^2)^ are the electron densities corresponding to the electron states, localized 
on one nucleus ("atom"), two nuclei ("molecule") and so on. As follows from fl33|) . fl35|) the 
total number of electrons Ne"^^ in the "localized" states can be represented as the sum of 
the total numbers of electrons in one-centrum ("atomic") N^[°'^\ two-centrum (molecular) 
N^!^^\ etc. states of the discrete spectrum. 

The ratios of the numbers N^'^"^\ ^f-rf* i depend on the thermodynamic parameters of the 
system under consideration. In the case, when there are only one-centrum ("atomic") states 
in the system under consideration, the total number of electrons in the "localized" states, 
according to ( I35|) .( l36l) is proportional to the number of nuclei: Ne°'^^ ~ N^. In this case, in 
the thermodynamic limit u'l"^^ ~ ric. 

The "delocalized" electron states describe the states in the continuous spectrum in the 
irregular field of a very large (infinite in the thermodynamic limit) number of nuclei. For 
the self-consistent description of the "localized" (in the one-centrum approximation) and 
"delocalized" electron states [20,21], electrons in the "delocalized" states interact not with 
"bare" nuclei, but with "ions" , which are nuclei dressed by "localized" electrons. 

In this case, to describe the " delocalized" states, it is natural to apply the perturbation 
theory to the interaction potential between "delocalized" electrons and "ions" [20,21]. Then 
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for the inhomogeneous density nh (r, {R-i}) of the " delocahzed" states, we can write 

^{deioc)^^^ {R=}) = nf''^'"'^ + 5nf''^"''\r, R^), (37) 

where 5?T,e'^^'°'^''(r, {R^}) is completely defined by the interaction potential between 
" delocalized" electrons and "ions" , and by the response functions of the homogeneous 
electron liquid, which is characterized by the temperature T and average density n^f^'"''^^ in 
the "delocalized" states. The examples of the use of the response functions in the density 
functional theory are given in [1, 4, 23-25]. 

Thus, we arrive at the conclusion that the problem of establishing of the explicit form of 
the density functional (fT8|) in the thermodynamic limit for the external potential of nuclei is 
the necessary step for finding the free energy functional F^. Even an approximate solution 
of this problem is very complicated and requires the use of the theory of disordered Coulomb 
systems. 



VI. CONCLUSIONS 



According to the above consideration we can assert the following: 

(A) The multi-component Coulomb system in the adiabatic approximation should be 
considered in the canonical ensemble. 

(B) The Hohenberg-Kohn theorem is valid for the free energy of inhomogeneous electron 
gas in an external field for non-zero temperatures. 

(C) The universality of the density functional for the free energy is meaningful with 
respect to the thermodynamic parameters of the system at the fixed external field and 
interaction potential. 

(D) In the DFT, there is a principal difference between the quantum-mechanical consid- 
eration of a finite number of particles in the infinite volume [Ne/V — 0, F — ?■ oo) and the 
quantum-statistical consideration in the canonical ensemble with the subsequent passage to 
the limit Ne oo, V ^ oo and n^e^ = N^/V — > constant. 

(E) Establishing the explicit form of the density functional for the external potential of 
nuclei in the thermodynamic limit is very complicated and requires of the use of the theory 
of disordered Coulomb systems. Descriptions of the "localized" and "delocalized" electron 
states are principally different. 
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